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The experimental STM images for the CDW phase of the blue bronze Rbo.sMoOs have been successfully 
explained on the basis of first-principles DFT calculations. Although the density of states near the Fermi level 
strongly concentrates in two of the three types of Mo atoms (Mon and Mom), the STM measurement mostly 
probes the contribution of the uppermost O atoms of the surface, associated with the MoiOe octahedra. In 
addition, it is found that the surface concentration of Rb atoms plays a key role in determining the surface 
nesting vector and hence the periodicity of the CDW modulation. Significant experimental inhomogeneities of 
the 6* surface component of the wavevector of the modulation, probed by STM, are reported. The calculated 
changes in the surface nesting vector are consistent with the observed experimental inhomogeneities. 



PACS numbers: 71.45.Lr, 68.37.Ef, 73.20.-r, 71.20.-b 

Low-dimensional molybdenum and tungsten oxides and 
bronzes have been the focus of much attention because of the 
charge density wave (CDW) and associated phenomena they 
exhibit 1 1]. The blue bronzes, A0.3M0O3 (A = K, Rb, Tl), are 
quasi-one-dimensional metals exhibiting a metal to insulator 
transition and are among the most intensely studied of these 
materials. Their crystal structure is built from M0O3 layers in 
between which the cations reside (Fig. [l}|2]. Despite many 
attempts, observation of the CDW in these materials by scan- 
ning tunneling microscopy (STM) has been elusive. Only very 
recently, high resolution STM images of an in- situ cleaved 
(201) surface of the rubidium blue bronze, Rbo.sMoOa, have 
been obtained at low-temperature and in ultra-high vacuum 
(UHV)(1. Both the molecular lattice and the CDW superlat- 
tice were observed simultaneously at temperatures well below 
the CDW transition temperature. 

Comparison of these images with previous first-principles 
density functional theory (DFT) calculations for the bulk0] as 
well as with experimental information of the bulk structure of 
the modulated phase] 5] is quite puzzling. For instance, some 
of the most intense features of the STM images are associated 
with the MoiOe octahedra (Fig. [l]), which are only weakly 
involved in the CDW transition according to the superlattice 
structural study 0]. Parenthetically, the Moi orbitals have a 
minor contribution to the states near the Fermi level accord- 
ing to the first-principles calculations] 4]. These and related 
observations prompted the present work. Here we report a 
first-principles study of the STM images of modulated and 
non modulated Rbo.sMoOs with special emphasis on the in- 
fluence of the alkali atoms at the surface. This is a key issue 
when trying to directly observe the CDW modulations in ma- 
terials like low-dimensional bronzes because it may affect the 
surface nesting vector and hence, the nature of the modulation 
observed. We also report new experimental results concern- 
ing inhomogeneities in the surface modulation wave vector 
providing support for the analysis. Finally, a clearcut under- 
standing of the STM images for the blue bronze emerges from 
this work. 
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FIG. 1: (Color online) (a) Idealized surface structure of Rbo.sMoOs 
in the (201) plane, (b) Idealized side view of the room temper- 
ature Rbo.sMoOs structure projected onto the plane perpendicular 
to the b axis. Closed circles are Rb atoms at the uppermost posi- 
tions of the "surface" (labeled 1) and empty circles are Rb atoms 
1.2 A below (labeled 2). The three highest octahedra with respect 
to the "surface" are the dashed squares indicated by the arrows. 
Their centers lie at levels 1.8, 2.4 and 3.5 A below the "surface". 

The present calculations were carried out using a numerical 
atomic orbitals DFT|6] approach, which has been developed 
and designed for efficient calculations in large systems and 
implemented in the SIESTA code|7]. We have used the gen- 
eralized gradient approximation to DFT and, in particular, the 
functional of Perdew, Burke and Ernzerhof|8]. Only the va- 
lence electrons are considered in the calculation, with the core 
being replaced by non-local norm-conserving scalar relativis- 
tic pseudopotentials|9] factorized in the Kleinman-By lander 
formyjD]. Nonlinear partial-core corrections to describe the 
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exchange and correlations in the core region were included 
for Mo| 11]. We have used a single-^ basis set including po- 
larization orbitals for Mo atoms, as obtained with an energy 
shift of 0.02 Ry|7]. We verified that the description of bulk 
bands for this and related bronzes using this basis size, espe- 
cially at Fermi level, is essentially the same as when using a 
split- valence double- C basis set including polarization for all 
atoms. The energy cutoff of the real space integration mesh 
was 300 Ry. Calculations for slabs of different thickness (con- 
taining from one to four octahedral layers) were carried out. 
For the superlattice we used the structure of reference 5 as- 
suming a commensurate value (0.75) of the b* component of 
the modulation vector. The Brillouin zone (BZ) was sampled 
using grids of (2x21x1) and (2x21x6) /c-points for the slabs 
and the bulk, respectively l[l2ll . The energy cutoff and /c-points 
values were tested against well converged values. 
The band structure for bulk Rbo.sMoOs contains two par- 
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FIG. 2: Band structure for: (a) bulk Rbo.sMoOs; (b) a slab 
preserving the bulk stoichiometry at the surface, and (c) a slab 
with a defect of Rb atoms (one every three) at the surface. In 

(a) r = (0,0,0), X' = (1,0,0) and = (0, |, 0) in 
units of the a'*, 6'* and c* reciprocal lattice vectors 0]. In 

(b) and (c) r = (0,0), X = (^,0) and Y = (0, |) in 
units of the corresponding rectangular reciprocal lattice vectors, 
tially filled bands (Fig. |2(a)| ). The CDW in this material is due 
to the interband nesting among these quasi-one-dimensional 



bands so that the CDW wave vector, qcDW^ is given by 
fefi + where k fi is the Fermi wave vector of band i 
lO [l4l lol . Since there are three electrons per unit cell to 
fill these bands, a Qcdw component along the chain direc- 
tion of 0.756* is predicted, which is the observed value at low 
temperaturel 16]. In order to appropriately model the (201) 
blue bronze surface we carried out calculations for slabs in- 
cluding different number of octahedral (and rubidium) layers 
as well as different distributions and concentration of surface 
rubidium atoms. There are three rubidium atoms per repeat 
unit of a layer, two of them (type 2, empty circles in Fig. |l(b)| ) 
are very near the octahedral layers and the third one (type 
1, full circles in Fig. |l(b)| ) is equidistant of the two layers. 
Among these interlay er Rb atoms, only the type 2 Rb atoms 
closest to the surface and the type 1 Rb atoms may remain at 
the surface after cleaving the sample. Since these atoms are 
expected to relax from their bulk crystallographic position, we 
optimized their position with respect to the surface. These po- 
sitions were the basis for all remaining calculations. 

The main conclusions of our study were: i) the number of 
octahedral layers used in the computations is irrelevant; ii) 
the key factor in controlling the shape of the surface bands 
near the Fermi level is the number of Rb atoms at the surface. 
Shown in Fig. |2(b)| is the band structure for a surface which 
preserves the stoichiometry of the bulk (i.e., 1.5 Rb atoms per 
repeat unit). Despite the folded-like shape of the bands, due to 
the fact that we used a unit cell twice larger along h in order 
to model the partial occupation of the Rb sites, the partially 
filled bands are nearly identical to those of the bulk. In con- 
trast, when the repeat unit of the cell at the surface is covered 
by just one Rb atom, i.e., 0.5 less than in the stoichiometric 
case, the band structure is noticeably different (Fig. |2(c)| ); the 
correspondig surface bands are shifted upwards with respect 
to those of the bulk. Calculations for the case of an excess of 
Rb atoms with respect to the stoichiometric situation led to an 
opposite band shift. After carrying out computations for sev- 
eral situations we conclude that different concentrations of Rb 
atoms at the surface generate a nearly rigid energy shift of the 
surface bands with respect to those of the bulk. This result has 
the important implication that the nesting vector at the surface 
changes with the Rb content. In fact, it is possible to infer the 
concentration of Rb atoms at the surface that produces a given 
surface nesting vector, see Fig.jSj 

Let us now search for the experimental consequences of 
this feature. As reported in previous work|3] we have ob- 
served with careful STM experiments the nearly commensu- 
rate value of the projection of qcDW onto the (201) plane. 
Hence, defining N=gr^* / 6* where gr^* is the h* component 
of qcDW^ N=0.75 is the bulk reported value and 1-N=0.25 
is the quantity reported by STM experiments (see part II and 
III in Ref. 3). However, it was mentionned that, on some 
optically flat terraces, 1-N was found to deviate from the 0.25 
value yielding inhomogeneities for gr^* . We report here exper- 
imental results concerning the inhomogeneities of q^^ mea- 
sured by STM in rubidium blue bronze. Three in- situ cleaved 
Rbo.sMoOs samples from the same batch were investigated 



Excess surface Rb atoms per repeat unit 

FIG. 3: Continuous line: calculated b* component of the surface 
nesting vector of the (201) surface of rubidium blue bronze versus 
the density of alkali atoms at the surface. The horizontal axis indi- 
cates the number of excess Rb atoms at the surface per unit cell, with 
the zero corresponding to the stochiometric Rb composition. Empty 
circles refer to the experimental values q* probed by STM (See text) 

with several mechanically sharpened Pt/Ir tips. All samples 
were prepared in a very similar manner. They were cleaved 
at room temperature in UHV and rapidly introduced into the 
cold STM head. All the STM measurements consisted of 
constant current mode topographical images and were per- 
formed at 77K or at 63K, well below the transition temper- 
ature (Tc=180K). Optically large and flat plateaus were care- 
fully selected to perform our measurements. Typical experi- 
mental conditions were ±450mV for the applied bias voltage 
and from 50 to 150pA for the tunneling current. Much care 
was taken in order to achieve molecular resolution and CDW 
resolution with scanning areas ranging from 20 x 20nm^ to 
50 X 50nm^. This was necessary to allow accurate measure- 
ments of the 1-N ratio for a typical 512 x 512 pixels image 
resolution. This 1-N ratio was directly extracted from the 2D 
Fourier transform of the STM image. At a given tip location 
several mesurements were always performed to ensure that the 
measured 1-N value was reproducible within an error bar of 
10%. 

It was found that optically distinct plateaux (of at least sev- 
eral 100 fim^ area) could yield distinct values of 1-N signif- 
icantly different from the 0.25 bulk value. Moreover, on the 
same plateau, different locations estimated to be at least sev- 
eral /im from each other yielded differences in 1-N values that 
were much greater than the typical error bar for a single loca- 
tion 1-N measurement, leading to clear inhomogeneities of 
the surface qcDW wave vector. On the contrary, displace- 
ments along 6 or a + 2c on the scale of tens of nanometers 
from a given position of measurement, did not lead to no- 
ticeable changes of gr^* . This shows that all measurements 
were performed far enough from CDW domain boundaries. 
On the same plateau the greatest change in 1-N value ranged 
from 0.21 to 0.32. As a result of our study the averaged 1- 
N values were found in the range 0.21 to 0.35, as indicated 
by the empty circles in Fig. [3l for the N value. According 
to the present calculation these inhomogeneities would cor- 
respond to an excess of surface alkali atoms ranging from 
about 0.1 to -0.35 per repeat unit. These predictions would 
then be consistent with the hypothesis that the distribution of 




FIG. 4: (Color online) (a) Constant current mode topographical 
image of 6.2 x 6.2 nm^ on (201) plane of Rbo.sMoOs at 63K 
(raw data image). The applied bias voltage is 420 mV and the 
set-up tunneling current is 110 pA. Molecular lattice and CDW 
superlattice coexist in the image. The three arrows indicate re- 
spectively the observed type I and II MoOe octahedra and the 
expected position of the MomOe octahedra. Associated profile 
along MoiOe octahedra idicated by left arrow (from Ref. f^). 
(b) (Color online) Calculated image and associated profile along 
the MoiOe octahedra for the modulated phase of Rbo.sMoOs. 



type 1 alkali atoms is responsible for the experimentally re- 
ported inhomogeneities. These inhomogeneities are consis- 
tent with photoemission and grazing incidence x-ray diffrac- 
tion results QX J8J, which both reported that about lOOK, 
N almost equals 0.75. This is because STM probes qcdw 
only at the uppermost layer of the compound and very locally 
at the nanometer scale inside a single CDW domain. On the 
other hand both x-ray and photoemission experiments probe 
deeper layers, which remain unaffected by the inhomogeities 
present in the first layer according to the present calculations, 
and over a much larger in-plane scale. This results in an aver- 
aging value of Qcdw which does not show any surface inho- 
mogeneities. 

We now turn to the analysis of the main features of the STM 
images for the modulated phase of Rbo.sMoOs in stoichio- 
metric conditions. As shown in Fig.|4l where the chains along 
b are readily visible, there is a very good agreement between 
the experimental and calculated images. The observed STM 
pattern inside the surface elementary unit cell consists of one 
well defined ball next to a more elongated continuous pattern 
along b and was attributedl 3] to the MojOe and MonOe octa- 
hedra respectively. In order to better understand these images 
we report in Fig. |5la plot of the iso-charge density integrated 
from the Fermi level to 0.5 eV above. Two features must be 
noted. First, there is essentially no contribution of the Rb 
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Oxygen of MOjO^ cations (Rb) 




FIG. 5: (Color online) Iso-Charge density (side view as 
in Fig. |l(b)t integrated from the Fermi level to 0.5 eV 
above. In the image in color the Rb, Mo and O atoms 
are shown as black, yellow and red balls, respectively. 

atoms. This provides computational support for the sugges- 
tion that the STM measurement is not sensitive to them| 19']. 
Second, the density is noticeable around the outer O atom of 
the MoiOe octahedra. These O atoms are the uppermost part 
of the surface, those of the MonOe octahedra staying approx- 
imately 0.6 A below. These two facts together easily explain 
why the brightest spots originate from the outer O atom of the 
MoiOe octahedra. In addition, the amplitude of the bulk verti- 
cal displacement of the outer O atom of the MojOe octahedra 
in the modulated structure used in the calculation is between 
three and four times smaller than the amplitude of the calcu- 
lated density profile along the MojOe octahedra (Fig. |4(b)| ). 
This means that the STM experiment is mostly measuring the 
differences in local density of states (LDOS) associated with 
these O atoms as a result of the existence of the CDW and not 
the differences in height of these atoms. The above mentioned 
apparent contradiction with previous works is solved by the 
results of Fig. [S] Around 42.5% of the charge density in this 
figure is associated to Mon, 22.1% to Mom, 1.9% to Moi and 
1.2% to the outer O of the MojOe. Consequently, the CDW 
modulation mostly affects the Mon and Mom atoms. How- 
ever, because of the non negligible participation of the outer 
O atoms of the MojOe octahedra, as a result of the strong 
hybridization between the Mo and O orbitals, the orbital mix- 
ing associated with the CDW modulation affects the LDOS of 
these atoms, leading to the differences in the profile of Fig. 
|4] According to the present results the more continuous path 
must originate from the outer O atoms of the MonOe octa- 
hedra but also from the Mon atoms which strongly partici- 
pate in the wave function. However the difference in height 
with respect to the uppermost part of the surface leads to the 



considerably less intense signal. Finally, the MomOe octa- 
hedra, which lie considerably deeper (i.e., approximatly 1.7 
A lower than the MojOe octahedra) are not visible at all in 
the STM image even though the CDW modulation strongly 
affects them. 

In summary, a combined theoretical-experimental approach 
has led to an in-depth understanding of the STM observation 
for the CDW in quasi-one-dimensional blue bronze. The puz- 
zling relationship between how CDW reveals in STM and 
x-ray experiments has been understood. The decisive role 
played by the surface Rb atoms, leading to experimental sur- 
face CDW wave vector inhomogeneities, has been clarified. 
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